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I.  INTRODUCTION 


Predicting  the  moment  exerted  by  a  liquid  payload  in  a  spinning  and 
coning  projectile  is  a  problem  of  considerable  interest  to  the  Army. 
Stewartson-1  considered  the  linear  problem  of  calculating  the  payload  moment 
through  the  use  of  separation  of  variables  and  eigenvalue  expansions  for*  an 
inviscid  liquid  in  a  cylindrical  cavity.  First  order  viscous  boundary  layer 
corrections  to  the  Stewartson  theory  were  carried  out  by  Wedemeyer2  and 
Murphy.3  A  method  for  calculating  the  linear  liquid  moment  using  the  full 
viscous  equations  with  boundary  layer  corrections  confined  only  to  the  endcaps 
was  also  presented  by  Gerber  and  Sedney.4  They  have  recently  extended  this 
theory  to  eliminate  the  boundary  layer  correction  at  the  endcaps.5 

Liquid  payloads  contained  in  a  highly  permeable  material  have  also  been 
of  interest  to  the  Army  for  some  time.  Laboratory  tests  and  flight  tests  have 
shown  that  a  highly  permeable  medium  can  significantly  reduce  the  spin-up  time 
of  a  liquid  payload. 6  Flight  stability  for  liquid  saturated  permeable 
payloads  has  also  been  examined  by  D'Amico.' 

This  report  extends  the  Stewartson  problem  by  considering  a  cylindrical 
cavity  filled  with  a  permeable  medium  that  is  impregnated  with  an  inviscid 
liquid.  A  further  modification  is  introduced  by  segmenting  the  cavity,  along 
the  symmetry  axis,  into  a  sequence  of  equal  length  cylinders.  Each  of  these 
cylinders  is  separated  by  impermeable  endcaps.  The  porous  media  is  modeled  by 
a  drag  term,  which  is  proportional  to  the  velocity  field,  added  to  the  linear¬ 
ized  Euler  equations.  This  analysis  examines  the  induced  liquid  moment  as  a 
function  of  parameters  found  by  Stewartson  plus  parameters  describing  the 
porous  media  and  the  number  of  segments  in  the  cylindrical  cavity. 

The  nomenclature  in  this  report  is  the  same  as  that  used  by  Murphy.3  In 
particular,  this  means  that  the  liquid  moment  can  be  represented  as  the 
complex  quantity: 


i4»c 

Transverse  Moment  =  mLaz<f>2  t  [C^m  +  i  <c  e 


(1.1) 


where 


mL  is  the  mass  of  liquid  in  a  fully-filled  cavity, 

a  is  the  maximum  radius  of  the  container, 

4>  is  the  inertial  spin  rate  of  the  container  along  the 
symmetry  axes, 

t  is  the  ratio  of  coning  rate  to  spin  <£c/$, 

Clsm  liquid  side  moment  coefficient, 

CLiM  is  the  liquid  in-plane  moment  coefficient, 

Kr  is  sin  «  is  the  precession  angle. 

^  c 


II.  LIQUID  MOMENT 

Two  coordinate  systems,  each  with  X-axis  along  the  projectile  symmetry 
axis,  are  used:  the  missile-fixed  (X,  Y,  Z)  system  and  the  non-rolling  XYZ 
system  with  the  Z-axis  initially  pointing  downward.  Introduce  an  earth-fixed 
axes  (Xe,  Yg,  Zg)  with  Xe-axis  in  the  direction  of  the  velocity  vector  and  Ze 
downward.  Let  a  unit  vector  in  the  direction  of  the  X-axis  have  earth-fixed 
components  (nXE,  nY£,  n^).  The  angle  of  attack  <*  in  the  (X,  Y,  Z)  system  is 

the  projection  on  the  XZ-plane  of  the  angle  between  the  X-axis  and  the  velo¬ 
city  vector.  The  angle  of  sideslip  3  is  the  projection  of  the  same  angle  onto 
the  XY-plane. 

The  kinematic  behavior  of  the  spinning  projectile  is  the  sum  of  two 
coning  motions:3 


«  =  6  +  1  a  "  Ki  e  1  +  K2  e 


(2.1) 


where 


ln  (VKjo)  -  ej  V 

♦j  =  ‘•’jO  +  Tj  * 


4>  =  <J>t 

(|  >0)  is  the  spin  rate  and  Kjn,  4. j Q  are  constants. 

III.  EQUATIONS  OF  LIQUID  MOTION  IN  POROUS  MEDIA 

The  following  analysis  models  the  steady  state  response  of  a  liquid 
flowing  in  a  porous  medium  confined  to  a  cylindrical  cavity  contained  in  a 
spinning  and  coning  projectile.  The  objective  of  this  theory  is  to  predict 
the  liquid  moment  resulting  from  coning  or  spiral  motion  which  is  specified 
by:3 


where 


5  =  K  e" 


-  ^iO 

K.Kj0e  ^ 


(3.1) 


j  =  1  or  2 


(Ej  *  ’>Tj 
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The  conservation  equations  governing  the  motion  of  the  confined  liquid 
are  the  continuity  equation  and  a  modified  Euler-momentum  equation.  The 
modification  is  one  that  is  commonly  used  to  describe  the  flow  of  liquid 
through  a  porous  medium.  This  consists  of  an  additional  term  given  by: 


where 


°r  *  -  7  <V  *  - 

d<p 


u  is  the  dynamic  viscosity 


(3.2) 


<  is  the  porosity  (dimensions  of  length2) 

4 

Vp  is  the  velocity  of  liquid  relative  to  the  porous  medium 
is  liquid  density 

4 

np  is  a  pressure  gradient  induced  by  resistance  of  the 

porous  media  to  fluid  flow.8  C  =  — -  —  is  a 

plK(? 

dimensionless  coefficient  which  is  a  measure  of  this 
pressure  gradient. 


Physical  reasons  for  using  such  equations  can  be  found  from  the  arguments  used 
in  establishing  Darcy's  Law.9 

Let  (r,0,x)  be  cylindrical  polar  coordinates  fixed  to  the  earth  frame  and 
(V,W,U)  be  the  corresponding  components  of  velocity.  For  small  angles,  the 
position  vector  of  any  point  in  the  projectile  has  components3 


r  =  r  -  K.  x  cos  (<p.  -  0) 

J  J 


x  +  K.  r  cos  ($. 

J  J 


0  +  0  (K?) 

J 


(3.3) 


(3.4) 


(3.5) 


where  tilda  (~)  quantities  are  measured  in  the  non-rolling  system.  Equations 
(3.3  -  3.5)  lead  to  the  velocity  components 


Vx  =  R{$(s  -  i)rK  es*  '  i6} 


(3.6) 


/p  =  _r{J(s  -  i )xK  es*  ‘  l6} 


(3.7) 


where 


Vfl  =  4>r  +  R { i 4>  ( s  -  i  )xK  e59  “  } 


R { }  =  [{}  +  {}  ]/2  =  Real  part  of  {} 


(3.8) 


Let  the  velocity,  V  ,  and  pressure,  p,  fields  have  the  form: 


V  =  q  +  4>  ex  xr 


P  =  p  + 


PL4>2r2 


(3.9) 


where  e  ,  efl,  e  are  the  unit  vectors  in  the  (r,e,x)  directions.  The 
r  o  x_^ 

variables  q  =  (v,w,u),  and  p  are  small  perturbations  of  0(Kj). 

Non-dimensionalizing  all  lengths  by  the  cylinder  radius  a  and  velocities 
by  a$,  plus  assuming  periodic  disturbances  of  the  form 

(v,w,u,p)  =  (v  (r,x),  w  (r,x) ,  p  (r.x))  es4”i9  (3.10) 


allows  the  continuity  and  momentum  equations  to  be  written  as 


1  .  1*  +  r  Hi  =  n 

dr  r  r  3x 


(s  -  i)  q  +  2  ex  x  q  =  -Vp  -  Cp[q  +  (s  -  i ) (x,-ix,-r) ]. 


The  components  of  Eq.  (3.12)  take  the  form 


(3.11) 


3.12) 


2w  +  -£2.  =  -  (s-i )  x  C 
dr 


2v  -  —  =  i  (s-i)  x  Cp 


tu  +  iH.  =  (s-i)  r  C 
3x 


(3.13) 


(3.14) 


(3.15) 


where  0<r<l,  h+f<x<h-f 


f  =  fineness  ratio 


Y  =  s-i  +  C, 


and  h  is  the  center  of  mass  location  along  the  symmetry  axis.  The  non- 
homogeneous  terms  of  Equations  (3.13  -  3.15)  suggest  using  the  transformations 


v  =  vH  - 


w  =  wh  + 


(s-i)  x  Cr 
(y  +  2i) 

i  (s-i)  x  Cr 
(Y  +  2i ) 


(3.16) 


(3.17) 


u  =  uH  + 


(s-i)  r  Cp 


(3.18) 


in  the  equations  of  motion  which  reduces  them  to  a  single  equation  for  the 
pressure  p 


32p  +  1  _  P_  _  _  (y2+4)  32p 


3r2  r  dr  r* 


Y2  3x< 


(3.19) 


Similarly,  the  solution  p  is  related  to  the  following  physical  quantities  of 
interest: 


r2  +  4 


[  zip  T  ie.  -] 


»„  -  -ji-  c  at ♦  2  ia  j 

Y2  +  4  r  3r 


u  =  .11 

H 

Y  3x 


(3.20) 


(3.21) 


(3.22) 


The  cylindrical  cavity  is  assumed  to  consist  of  N  chambers  with  circular 
cross  section  and  height  a  =  2f/N,  each  separated  by  impenetrable  endcaps  as 
depicted  in  Figure  1. 


The  boundary  conditions  for  Eq.  (3.19)  are  the  normal  velocity  of  the 
fluid  at  the  boundary  must  equal  the  normal  velocity  of  the  wall  and  all  flow 
variables  remain  finite.  This  means: 


i£  -  o 

3x 


atx  =  h-  f  +  n  ( — ) 
N 

n  =  0,1,... N 


(3.23) 


Y  l£  =  .  2xs  (s-i)  (s-3i+C  ) 
dr 


at  r  =  1 


(3.24) 


p  =  0  at  r  =  0 


(3.25) 


where 


=  p  +  xr  (s-i )2. 


(3.26) 


This  boundary  value  problem  is  similar  to  the  problem  considered  by 
Stewartson1  and  can  be  shown  to  have  the  solution 


p  =  -  xr  (s-i)2  +  xs  (s-i)  [2 ( h-f )  +  (2n+l)  A]  r  + 


16fs  (s-i )  (Y-2i ) 

N  k^d 


Jj(Xkr)  cos  ij-  (x-h+f-nA) 
k2  [(21+y)  (Xk)  -  r  xk  JQ  (Xk)] 


(3.27) 


where  xj^  = 


(y2+4)  N2  k2  tt 2  . 
Y2  4  f 2 


N  =  1,2,3,... 
n  =  0,1 . N 


and  the  Jn  are  the  Bessel  functions  of  the  first  kind.  If  Y  is  pure 

imaginary,  the  transcendental  denominator  in  Eq.  (3.27)  will  have  zeroes. 

This  resonant  condition  occurs  whenever  t  e  +  C  =  0  and  i  is  one  of  the 
inviscid  Stewartson1  eigenf requencies. 


IV.  LIQUID  MOMENT  EQUATIONS 

The  moment  induced  by  the  liquid  contained  in  the  segmented  cavity  is 
calculated  from  the  time  derivative  of  the  angular  momentum  field.  Non- 
dimensional  izing  this  moment  with  2-jt a 3  f  allows  the  moment  in  the  YX  plane 

to  be  expressed  as  a  single  complex  quantity3 


My  +  i  Mz  =  T  cLM  K  eJT 


(4.1) 


where 


^LM  =  ^LSM  (T,e»^r»^)  +  ^  ^LIM  (T*e»^r»^)* 


The  unit  vectors  for  the  earth-fixed  cylindrical  coordinates,  (e-  ea,  e  ), 

r  ti  x 


can  be  written  in  the  same  complex  notations  in  terms  of  (e^,  ez,  ex) 


e  =  e  cos  6  +  e  sin  6  <=>  e 
r  y  z 


*  *  *  j  Q 

eA  =  -  e,  sin  6  +  e,  cos  e  <=>  ie  . 
a  y  z 


(4.2) 


(4.3) 


Substituting  these  into  the  moment  integral  and  using  the  Reynolds  Transport 
Theorem10  gives  the  following  expression  for  the  liquid  moment  coefficient: 

Nl-  rh-f+  (n+1  )A  ,1 

T  CLM  =  ~  XI  |  «*x  +  ?er^  x  [(S  "  +  2ex  x  ^  (4*4) 

2^  n=0  h-f+nA  •'0  -r~?  ~  , ~ 

-  i [r2  -  2x2]}  r  dr  dx. 

Using  Eqs.  (3.16-3.18  and  3.20-3.22)  with  the  aid  of  Eq.  (3.27)  permits 
writing  the  last  equation  as 

i  f4[l-C  (s-i)]  [(sM)  (3h2+f  2)]Cp  3[1+C  (s-i)](s-i)2Cr 


1 2h 2  -  4f 2 


+  3  +  .  -211—  (111 

(y+2i) 


1)  (3h2+f 2 


128  i  (s-i)  (y-2i )s2f2 


^+2i)N4lt4  k~d  k4K^2i)  Jx(Xk)  -  YXk  00(Xk)] 


N  =  1,2,3,... 


Murphy3  derived  the  frozen  liquid  values  of  C,  cM  and  Ci  TM  and  these  are  given 


'.4.1 


=  "  l 
2 


1  -  T  1  + 


4  (3  h2 


r1]] 


;  =  1  +  T.(e2  U  [3  +  12h2  +  4f2]  . 

2  12 


(4.6) 


(4.7) 


For  very  slow  motion  (t  ♦  0)  the  liquid  moment  should  approach  these  values. 
Equation  (4.6)  shows  that  the  limiting  value  of  CLIM  as  t  0  equals  (e  + 

i)/2.  This  agrees  with  the  forzen  liquid  results  when  t  =  0.  When  Cr  ♦  "the 

liquid  should  act  like  a  frozen  liquid  for  all  values  of  t.  It  is  easy  to  see 
that  the  liquid  moment  coefficients  given  by  Equation  (4.6)  approach  (in  the 
limit  of  Cp  +  «  for  Ak  not  a  resonance  value)  the  values  for  a  frozen  liquid. 

The  same  limits.  Equations  (4.6  -  4.7),  are  also  found  when  Cr  and  t  are 
arbitrary  hut  N  becomes  infinitely  large. 

V.  SOLUTION  METHOD 

The  equations  of  sections  III  and  IV  need  to  be  solved  over  a  wide  range 
of  t  for  specific  values  of  0  ,  f,  and  N.  This  requires  calculating  the  ratio 

J0(Ak)AMV  the  zero  and  ^rst  order  Bessel  functions.  For  values  of  k 

such  -.hat  |Ak|  <  50,  the  above  ratio  is  obtained  by  simply  dividing  the  values 

found  from  power  series  expansions  of  each  Bessel  function.  For  larger  values 
of  |Ak|,  an  asymptotic  expansion  of  the  Bessel  function  ratio  was  used.11 

Equation  (4.5)  is  then  used  to  find  values  of  C|_<^  and  Experience  has 

shown  that  k  <  20  i s  sufficient  to  produce  converged  solutions.  All 
computations  were  carried  out  on  a  VAX-8600  computer. 


VI.  DISCUSSION 

Figures  2  and  3  present  plots  of  C|_ 5^  and  C[_j^  as  functions  of  frequency 

for  f  =  1.5,  Cr  =  3.0,  N  =  1,  and  e  =  0  and  0.02.  For  zero  damping,  the 

maximum  side  moment  is  slightly  larger  than  the  maximum  side  moment  for  a 
small  amount  of  undamping  given  by  e  =  0.02.  The  in-plane  moment  remains 
relatively  unchanged  due  to  the  presence  of  the  same  undamping.  This  indi¬ 
cates  that  CLM  is  insensitive  to  e  for  e  near  zero. 

It  was  mentioned  earlier  that  C|_jM  and  CLjM  approach  the  limiting  values 

for  a  frozen  liquid,  Eqs.  (4. 6-4. 7),  as  Cr  *  0.  An  example  of  this  for  f  =  2, 

N  =  1  and  e  =  0  is  exhibited  in  Figures  4  and  5.  Similar  results  are  shown  in 
Figures  6  and  7  for  an  increasing  number  of  chambers  N  and  a  fixed  Cp  =  3. 


Generally,  it  has  been  found  that  the  frozen  liquid  limit  is  approached  quite 
rapidly  with  increasing  values  of  N  for  typical  values  of  Cp. 


If  tq  is  a  Stewartson  eigenfrequency  for  the  particular  values  fQ,  kQ, 
and  Nq,  then  a  resonant  condition  will  occur  at  t  =  t0  whenever 


k q ,  k  =  1,3 ,5, ■  • » 


N0,  N  =  1,2,3,... 


(6.1) 


For  an  example  of  this  phenomenon,  consider  the  case  of  kQ  =  1,  f  =  2,  NQ  = 
1,  Cr  =  3,  and  e  =  0.  These  parameter  values  make  the  ratio  of  Equation  (6.1) 
equal  to  2  and  the  first  eigenfrequency  tq  =  0.5102.  Figure  8  shows  CL$M  as  a 
function  of  t  for  k  =  1,  f  =  10,  and  N  =  1,5,10.  The  increase  in  CL5M  for  N  = 
5  is  due  to  the  eigenfrequency  tq  =  0.5102  since  the  ratio  in  Eq.  (6.1)  is 
again  equal  to  2.  Hence  for  a  fixed  t,  an  increase  in  N  can  cause  CL^M  to 
increase.  The  eigenvalues  for  the  N  =  1,10  cases  lie  outside  the  range  of  t 
given  in  Figure  8.  CLSM  therefore  decreases  monotonically  to  the  frozen 
liquid  value,  which  is  zero  for  e  =  0,  when  N  becomes  sufficiently  large. 

The  155mm  M825  projectile  is  shown  in  Figure  9.  This  projectile  carries 
a  canister  that  is  loaded  with  white  phosphorous  (WP)  impregnated  felt  wedges. 
Four  angular  ribs  produce  longitudinal  quadrants  within  which  felt  wedges  (116 
per  canister)  are  loaded.  Aluminum  foil  spacers  are  located  between  each 
wedge  (as  shown  in  Figure  9)  and  could  produce  a  compartmentalization  or 
segmentation  as  modeled  in  this  theory.  When  the  WP  is  liquid  (temperatures 
above  44  deg  C),  flight  instabilities  have  been  recorded.  Recent  flight  tests 
have  shown  that  the  flight  stability  of  the  projectile  can  be  improved  by 
using  felt  wedges  with  outer  diameters  that  produce  interference  fits  when 
loaded  into  the  canister.12  The  model  given  in  this  report  is  directed  at  a 
fundamental  understanding  of  the  payloads  used  in  the  M825.  The  inclusion  of 
a  drag  force  on  the  liquid  caused  by  the  permeable  media  provides  a  first  step 
in  understanding  the  physics  of  these  payloads. 

Scheidegger13  tabulated  permeabilities  of  various  substances  and  stated  a 
range  for  hair  felt  as  8.3  x  10”6  cm2  <  x  1.2  x  10"5  cm2.  For  an  M825-type 
payload,  nominal  values  for  spin  rate  and  kinematic  viscosity  are  100  Hz  and 
0.015  cm2/sec.  Hence,  a  median  drag  coefficient,  Cr,  is  2u/15  or  approximate¬ 
ly  1/2.  Experiments  should  be  conducted  to  determine  the  permeability,  k,  for 
the  type  of  felt  and  packing  used  in  the  M825  projectile.  Investigations  must 
be  conducted  to  verify  if  a  linear,  homogeneous  formulation,  using  Darcy's 
Law,  is  appropriate  for  spinning  liquids.  The  present  analysis  extends  prior 
theories  for  liquid  payloads  and  produces  a  simple  model  from  which  a  clear 
physical  understanding  is  obtained.  Additions  to  the  present  effort  could 
include  the  effects  of  fluid  interaction  with  solid  boundaries,  radial 
variations  in  Cr  and  the  nonisotropic  character  of  Cp. 


VII.  CONCLUSIONS 


The  liquid  moment  coefficients  are  computed  for  cylindrical  cavities 
which  are  fully-filled  with  a  permeable  medium  and  impregnated  with  an  invis- 
cid  liquid.  These  coefficients  reflect  the  results  of  segmenting  a  given 
cavity  into  N  chambers  with  uniform  height.  The  coefficients  are  then  calcu¬ 
lated  as  functions  of  coning  frequency,  fineness  ratio,  and  the  parameter  Cp 

used  to  represent  the  physics  of  a  liquid  flowing  through  the  permeable 
medi urn. 

If  the  range  of  coning  frequencies  does  not  contain  any  eigenfrequencies, 
then  the  calculations  plus  theory  indicate  that  the  side  moment  approaches  the 
value  for  a  frozen  liquid  as  Cr  +  0  and/or  N  becomes  large.  However,  a  signi¬ 
ficant  increase  in  the  side  moment  can  occur  for  a  particular  N  by  choosing 
the  fineness  ratio  such  that  an  eigenfrequency  moves  into  the  frequency  range 
of  interest. 
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LIST  OF  SYMBOLS 


a 

CLI«j 


Maximum  radial  distance  of  the  liquid-filled  container  (for  a 
cylinder,  the  radius,  for  a  spheroid,  the  radial  semi-axis) 

Imaginary  part  of  C|_M_;  coefficient  representing  the  liquid 
moment  that  causes  rotation  in  the  plane  of  exp  i$),  j  =  1,2 

J 


1 


Clsmj 

Cr 

°r 


Real  part  of  C|_M.;  coefficient  representing  the  liquid  moment 
that  causes  rotation  out  of  the  plane  of  exp  ( i ^ ) ,  j  =  1,2 

J 

Drag  coefficient 

Drag  term  due  porous  media 


Unit  vectors  along  the  earth-fixed  cylindrical  axes 


Unit  vectors  along  the  earth-fixed  Cartesian  axes 


f  Fineness  ratio,  c/a,  for  a  cylinder 

h  Distance  of  center  of  gravity  along  symmetry  axis  from 

geometric  center 


k 

K 


Kjo 


mly  +  mlz 


Axial  wave  number 

Kj0e14>j0  (j  =  1,2) 
e-j  Ti  ‘l’ 

Kjoe  J  J  (j  =  1,2) 

Value  of  Kj  at  t  =  0 

Mass  of  the  liquid  in  the  container  2na2cp,  for  a  cylinder; 
4na2cp^/3  for  a  spheroid) 

Transverse  liquid  moment  in  the  aeroballistic  nonrolling  system 


N 

nXE,nYF’nZE 

P 

r 

RO 


Number  of  products  in  the  X  factor,  Eq.  (3.19) 

Earth-fixed  components  of  a  unit  vector  along  the  x-axis 
Non-dimensional  pressure  perturbation  as  a  function  of  r  and  x 
Radial  coordinate  in  the  earth-fixed  cylindrical  (x,r,e)  system 
Real  part  of  {} 
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WVe 


X  Y  7 
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uj  *  (),j 


j  =  1  or  2 


Non-dimensional  velocity  perturbation  components  in  the  earth- 
fixed  cylindrical  (x,r,e)  system 

Non-dimensional  liquid  velocity  perturbation,  Eq.  (A. 3) 


Components  of  V  in  the  earth-fixed  cylindrical  (x,r,e)  system 


Axial  coordinate  in  the  earth-fixed  cylindrical  (x,r,9)  system 

Missile-fixed  axes,  the  X-axis  along  the  projectile's  axis 
of  symmetry 

Aeroballistic  non-rolling  axes,  the  Z-axis  initially  downward 

Earth-fixed  axes,  the  Xe-axis  along  the  velocity  vector, 

Ze  downward 


(K./K.)/$->  non-dimensional ized  damping;  j  =  1,2 
J  J  J 


Azimuthal  coordinate  in  the  earth-fixed  cylindrical 
(x,r,e)  system 

Kinematic  viscosity  of  the  liquid 


Defined  in  Eq.  (3.27) 

Liquid  density 

$./<£>  non-dimensional  ized  frequency;  j  =  1,2 
J 

Spin  rate  with  respect  to  inertial  axis,  assumed  positive 
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